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ABSTRACT 

We derive the rest-frame i^-band luminosity function for galaxies in 32 clus- 
ters at 0.6 < z < 1.3 using deep 3.6/zm and 4.5/im imaging from the Spitzer 
Space Telescope InfraRed Array Camera (IRAC). The luminosity functions ap- 
proximate the stellar mass function of the cluster galaxies. Their dependence on 
redshift indicates that massive cluster galaxies (to the characteristic luminosity 
M^) are fully assembled at least at z ~ 1.3 and that little significant accretion 
takes place at later times. The existence of massive, highly evolved galaxies at 
these epochs is likely to represent a significant challenge to theories of hierarchi- 
cal structure formation where such objects are formed by the late accretion of 
spheroidal systems at 2; < 1. 

Subject headings: galaxies: luminosity function, mass function - galaxies: forma- 
tion and evolution 

1. Introduction 

Clusters of galaxies are important for studies of galaxy formation and evolution, because 
they contain a volume limited population of galaxies observed at the same cosmic epoch. 
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They therefore provide a well-defined sample of objects to cosmologically significant lookback 
times whose member galaxies can be identified by simple counting statistics, without need 
for extensive spectroscopic surveys or multi-color data. 

One important characteristic of early type galaxies in clusters is that they are known 
to follow tight color-magnitude relations, which appear to be universal and to have very 



small intrinsic scatter to the highest redshifts yet observed ( Visvanathan fc Sandage 



Bower et al. 


1992; 


Stanford et al. 


1995 


Holden et al. 


2004: 


Mei et all 


2006a 


b) 



2003 



Lopez-Cruz et al. 



1977 



2004 



bj). Together with the conventior ial interpretation o f 



the color- magnitude relation as a mass-metallicity correlation (e.g., iTrager et al.l |2000| ). 
this implies that the majority of the stellar populations in early-type cluster galaxies were 
formed via rapid dissipative starbursts at z > 2. Fundamental plane studies of high red- 
shift cluster galaxies also support this conclusion , at least for the rn ore massive objects 
(Ivan Dokkum fc Stanford! I2OO3I : IWuvts et all I2OO4I : Iflolden et aP boosh . al though the low 
mass galaxies seern to have undergone more ext ended star formation histories (jPoggianti et al. 
200l[ iNelan et allboosl : Ijorgensen eraPboOsh . 



Theoretically, the existence of such massive and old galaxies at high redshift should 
represent a severe challenge to models where galaxies are assem bled hierarchically, from a 
sequence of major mergers at progressively lower redshifts (e.g., IColesI l2005l : ISpringel et c 
2OO5I : lBaughll2006l for recent reviews). It is not possible, however, to exclude by spectropho- 



tometry alone, that these galaxies are assembled from sub-units whose star formation has 
already c eased, but which are not accreted until l ater times (similar to the so-called 'dry' 



van 



Dokkum et aDll999l : iTran et al.l[2005l ) . This is assumed to be the main channel 



mergers 

by which spheroids grow at 2; < 1 in the hierarchical scenario. 

On the other hand, if galaxies are formed via mergers, we should observe a steady 
decrease of the mean stellar mass in galaxies as we go to higher lookba ck times and the most 



massive members of th e merger tree branch into ever smaller twigs (jPe Lucia et al.l 12006 



Maulbetsch et al.ll2006l ). While it is generally difficult to measure galaxy masses, the K- 



band luminosity function is believed to provide an adequate surrogate (IKauffmann fc Chariot 



19981 ). as the rest frame H 01 K luminosity of galaxies is seen to correlate well with s tellar 
and even dynamical mass for l ocal samples ( Gavazzi et al. 19961: Bell &: de JongI 2001) ari d 



even 



for high redshift galaxies ( Drory et al.lboO^ : IPapovich et al. 



2OO5I : ICaputi et al.ll2006h . 



In our previous work ( De Propris et al]|l999l ) we showed that the observed (ground- 
based) X-band luminosity of galaxies in clusters was consistent with pure passive evolution 
of objects formed at high redshift and argued that this implied that the majority of the 
stellar mass was completely assembled by at least z = 0.9. More recent luminosity function 
studies have essentially confirmed and extended this picture of early galaxy assembly in 
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clusters ( 


Kaiisawa et al. 2000: T 


^akata et al. 


2001: 


Massarotti et a 


. 2003 


Kodama & Bower 


2003 


; Toft, Soucail & Hiorth 2003; 


Ellis & JonesI 


20o4 


Toft et al. 


2004 




Bremer et al.l 


2006: 


Lin et al. 


2006: StrazzuUo et al. 


2006 


)■ 


Andreon ( 


2006 


recently derived a composite 3.6/im 



luminosity function for galaxies in clusters in the XMM-LSS survey (at a mean redshift of 
0.5), finding that the results are consistent with the previous ground-based results. 

Ideally, we would wish to carry out this experiment in the rest frame K-hand, as even 
the ground-based K-ha nd starts starts to sample the rest-frame op tical ai z > 1. T he 
Spitzer Space Telescope ( IWerner et al.ll2004j ) InfraRed Array Camera (IFazio et al.ll2004al ) is 
now capable of obtaining panoramic (~ 5' x 5') images at A > 3/xm with fxJy sensitivity 
and allows us to study the rest-frame i^'-band luminosity function of cluster galaxies at high 
redshift. 



Here, we present a study of 32 clusters up to z = 1.3 in both the 3.6;um and 4.5^m bands 
and derive the evolution of the rest-frame K-hand galaxy luminosity function, which is a 
close proxy for the stellar mass function. We adopt the cosmological parameters Qm = 0.3, 
l^A = 0.7 and Hq = 70 km s^^ Mpc^^ 



2. Data Reduction and Photometry 



The sample consists of 32 clusters at 0.6 < z < 1.3. Data were obtained with IRAC 
in all four filters, using 5 dithered frames of 200s each. Here we discuss observations in the 
3.6/zm and 4.5/xm filters, which map more closely to the rest-frame K-band at the redshifts 
of our clusters. 

Table 1 shows a list of clusters and some re levant properties. Most of the sample comes 
from the ROSAT Deep Cluster Survey (RDCS - iRosati et al.lll998l ). while a few others derive 
from other X-ray or optical surveys (see table for details), but the target sele ction is somewhat 
hetero gene ous, especially f o r the h igher redshift objects. On the other hand, iDe Propris et al. 
( 120031 ) and iPopesso et al.l (120051 ) have shown that the S-band galaxy luminosity function 
does not depen d on cluster properties s uch as the velocity dispersion, Bautz-Morgan type 
or richness and iDe Propris et al.l (119991 ) found no difference in the ground-based i^-band 
luminosity function of clusters selected by density and X-ray luminosity. 

The IRAC data were reduced following standard procedures. The raw (BCD) data were 
first corrected for known IRAC artifacts associated with bright stars (mux-bleed and column 
pulldown). Then the SSC (Spitzer Science Centre) MOPEX software was used to mosaic the 
individual frames into a registered mosaic, with cosmic rays removed. This mosaiced image 
for each cluster and for each band has a pixel scale which is 1.414 times smaller than the 
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original 1.2" IRAC pixel scale, and the orientation rotated by ~ 45^^ 



Photometry was carried out using the Source Extractor software (IBertin fc Arnouts 



19961 ). We experimented with various values for the background level and the deblending 
threshold, because Spitzer images have relatively poor angular resolution (1.7" FWHM for 
stellar sources) and because our fields are moderately crowded. 

In order to deal with the moderate crowding, we checked that the poorer resolution 
of Spitzer data does not significantly affect our detection and photometry. We verified 
the detections visually, both on the original image and on the aperture image produced 
by Sextractor. We also compared photometry in 2", 3" and Kron apertures, extrapolated 
to total magnitudes, to check that objects were properly deblended. These tests provide 
confidence that our photometry is not significantly affected by the crowding, although to 
fully address this issue will require higher resolution imaging. 



Number Counts for Cluster Galcixies 



We chose to measure magnitudes in fixed 3" (radius) apertures, which we re calibrated on 
to the Vega system and extrapolated to infinity following iFazio et al.l (l2004bf) to produce to - 
tal magnitudes. This is done for consistency with the apertures used by lFazio et al.l (l2004bl ) 
to derive galaxy number counts in IRAC bands, which we use for background subtraction. 
We then counted objects (stars and galaxies) in 0.5 magnitude wide bins within a circular 
aperture of radius 1 Mpc, centred on the brightest cluster galaxy (where the cluster density is 
higher with respect to background). These systems tend to lie at or near the peak in galaxy 
density and the centre of the X-ray isophotes. Ideally, we would wish to choose an aperture 
based on the cluster structural parameters (e.g., r2oo) but the IRAC field size is not suffi- 
ciently large to derive a reliable profile for the cluster galaxy distribution. Since it is unlikely 
that the few bright galaxies at large cluster-centric radii may affect the luminosity function 
parameters significantly, and since there is little evidence that the lu minosity function varies 



with distance from the clus ter centre, at least for bright galaxies (jPe Propris et al.l 12003 



Lin, Mohr & Stanford 2004. iPopesso et al.ll2005l ). our choice of a 1 Mpc aperture should not 
affect our conclusions. 

We estimated the contrib ution of background galaxies to the observed counts by using 
the 3.6/im and 4.5/im counts of iFazio et al.l (j2004bl ). We fitted a low-order polynomial to the 
literature counts to smooth the effects of large scale structure along the lines of sight of the 
background fields. Error s were assumed to b e Poi ssonian, while the c lustering contribution 
was calculated following iHuang et al.l (119971 ) and iDriver et al.l ( l2003l ) . The Poisson errors 
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for the cluster counts and the background contribution and the clustering errors for the field 
contribution were then added in quadrature as appropriate. 

Because of the low resolution of Spitzer data, we are not able to discriminate easily 
between stars and galaxies. There are no published star count models for Spitzer passbands. 
We estimated the stel lar contribution us ing the predicted L band counts from the Besangon 
m odel of the Galaxy (IRobin et al.ll2003l ). These give a good fit to the star counts reported 
bv IPazio et aD J2004b[ ) in the Extended Groth Strip and QSO 1700 fields. 



Table 2 shows the raw number counts, estimated background contribution, stellar con- 
tamination and corrected number of (statistical) cluster members to the limiting apparent 
magnitude we use for both IRAC bands. The limiting magnitude is designed to reach the 
same absolute magnitude in all clusters (in two broad redshift ranges - see discussion be- 
low), such that the cluster counts in the faintest bin are still significantly above the predicted 
contamination. At the same time, the brighter magnitude limit (typically around 18*'' mag- 
nitude) reduces the effect of crowding, which is most significant for the fainter galaxies. This 
limit is found to lie below the knee of the luminosity function and in the regime where the 
counts are fitted by a power law. 

There are some objects where the level of contamination from foreground stars is high 
or which have low number counts. These objects are identified in Table 2 and not used in our 
analysis. In practice, we choose objects where the residual cluster counts (after removal of 
background galaxies and stars) are higher than 50 in the 3.6;um band for the z < 0.9 sample 
(where this band is closer to the rest-frame K) and larger than 25 for the z > 1.1 sample in 
the 4.5//m band (i.e., where this passband better probes the rest-frame K) - C10848. 9+4452 
is an exception to this, as we have only 1/4 of the field of other clusters. We remark that 
there is no evidence for the actual existenc e of our h i ghest redshift target (QSO1215-00), 
suggesting that the structure identified by iLiu et al.l (120001 ) consists of a small group or 
filament. 

The actual counts for each cluster suffer from small number statistics. Rather than 
recovering the luminosity function from Bayesian methods (e.g., Andreon, Punzi & Grado 



2005) we use composite luminosity functions, in order to average errors out (c.f. [Andreon 
20061 ). We create composite luminosity functions for clusters in two redshift bins, a t 0.6 < 
z < .9 and 1.1 < 2; < 1.3, in both bands, following the procedure described by IColless 
(119891 ). We bin galaxies in 0.5 magnitudes wide absolute magnitude bins in rest-frame 



K, adopting the cosmolo gy specified above and a ^-correction derived using the models of 
Bruzual fc CharlotI (120031 ) to transform from the observed IRAC bands to rest-frame K. We 
choose to sar aple these two redshift bins for the following reasons. Most previous studies, 
starting from lDe Propris et al.l (119991 ) have studied clusters a.t z < 1; only recently have ade- 
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quate cluster samples at z > 1 become available (e.g.. lToft et al.ll2003l . 12004 IStrazzullo et al. 



20061 ). The two redshift bins we study sample the rest frame K-band luminosity function 
of galaxies in these two regimes, i.e. both the reasonably well studied z < 1 sample and 
the more recent clusters at higher redshift. Furthermore, the 3.6 ^m band maps closely to 
rest frame K for the 0.6 < z < 0.9 interval, while the 4.5 fim band does the same for the 
1.1 < 2; < 1.3 regime. 

The fc-correction used above assumes a solar-metallicity single stellar population formed 
at 2; = 3 and with star formation declining exponentially with an e-folding time r) of 0.1 
Gyr, and is computed independently for each Spitzer band, which is thus transformed to rest- 
frame K. This is done (rather than more complex approaches using both IRAC bands to 
derive the rest-frame K luminosity) to present the data more directly and with a minumum 
of model dependencies. The /c-corrections used are presented in Table 3 (the full version is 
available electronically, while we only show the first ten lines of the table for guidance here). 
We experimented with several 'reasonable' values of r from instantaneous star formation to 
1 Gyr e-folding time bursts and found that this makes little difference to the actual value of 
the /^-correction. 

The resulting composite luminosity functions are fitted with a Schechter function, using 
a downhill simplex algorithm. Figure 1 shows the data in each bin and the best fitting 
luminosity functions. Table 4 shows the derived values for the luminosity function in 
both bands. The errors in M* are marginal la errors derived by fixing the values of all other 
parameters at their 'best' value. The derived a is also shown in Table 4, but we caution that 
the fit to the faint-end slope is very uncertain. 



Discussion 



Figure 2 shows the rest-frame K band M* for our composite clusters, together with 
previous ground- based K band data, corrected to rest-frame K following the same procedure 
as above, and a iBruzual fc CharlotI (120031 ) model, with solar m etallicity (Salpeter Initia l 
Mass Function and Padova 1994 isochrones, as recommended by iBruzual fc CharlotI l2003l ) 
and variable formation epoch, with r = 0.1 Gyr. The actual of choice of r makes a difference 
only to the level of a few hundreds of magnitudes. Note that we are not attempting to actually 
fit models to the data shown in Fig. 2, but we are showing models with representative star 
formation histories to obtain upper limits to the epoch of mass assembly for these cluster 
galaxies. 



The results shown in Fig. 2 imply that the majority of the stellar mass in elliptical 
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galaxies is already assembled at least at z = 1.3: this is a strong upper limit to the epoch 
of galaxy formation in that the majority of the merger episodes (if any) must have taken 
place prior to this epoch. Taken at face value, these Spitzer data may indicate that the 
epoch of star formation in these objects is somewhat more recent (1.5 < z < 2) than 
indicated by some previous studies of t 



1998; Blakeslee et al. 2003; Holden et al. 



le color-magnitude rela ti on (iStanford et al. 



12004 



2004; Mei et al. 



but see 



2006a 



b; 



J0rgensen et al. 



Holden et al.l l2006l ) and 



1995 



20051 ). This may be due 



the fundamental plane (IWuyts et al. 
to the fact that, at ~ M*, the population includes a fraction of lenticular galaxies, whose 
star formation histories are more extended than the brighter elliptic al galaxies. Otherwise, 



our results are consistent with recent work (e.g. iHolden et al.ll2006l ) showing that massive 
ellipticals are in place at 2; ~ 1 but are more general, as we make no morphological selection. 



The latest hierarchical models (e.g., iDe Lucia et al.l l2006l ). including AGN feedback 



succeed in pushing back the epoch of major star formation to z > 2 for the most massive 
objects, but still require most of the actual mass assembly to take place through dry mergers 
at later epochs; 50% of the mass in more massive galaxies (M > 10^^ ^0) is assembled 
at z < 0.8, while the lower mass objects (M > 4 X 10^ Mq) may be formed at higher 
redshifts. This is not what we observe here, where the vast majority of the stellar mass 
in massive (approxir nately L* or greate r, equivalent to a mass of ~ 10^^'^ Mq, using the 



calibration shown in lGavazzi et al.lll996l . which corresponds more closely to the virial mass 
- using stellar masses this is approximately 10^^'^ Mq) galaxies appears to be in place at 
z = 1.3. 

One possible caveat is that the theoretical models refer to the average 'field' environment, 
while we are observing massive clusters where the main process of hierarchical merging and 
collapse may have taken place at earlier epochs, as they lie in overdense regions. However, 
Maulbetsch et al.l (120061 ) use a high-resolution N-body simulation to simulate how the mass 
assembly histories of galaxy-size haloes depend on environment and show that aX z = 1 the 
mass aggregation rate is 4 times higher than at present, and independent of environment, 
while galaxies in the densest (cluster-like) environments at z > 1 undergo more rapid mass 
accretion. This suggests that we should be witnessing a much stronger evolution of the mean 
galaxy mass than observed here, even for cluster environments. 

Our results are therefore largely inconsistent with the hierarchical picture. Not only 
are the stellar populations of these galaxies formed at high redshift (see Introduction), but 
they are also assembled into galaxies at comparatively high lookback times, arguing that 
star formation takes place in situ, in a manner reminiscent of the earlier monolithic collapse 
picture. Recently, it has b een shown that this behavior largely hold s for field ellipticals as 
well, at least to z ^ 0.65 (IRoseboom et al.ll2006l ; IWake et al.ll2006l ). Similarly, i^-selected 
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studi es in the field have also found a similar anti- hierarchical b ehavior (e.g., ICimatti et al. 



2004 from the K20 survey - see also the review by iRenzinil |2006| for the observed 'top down' 
buildup of massive ellipticals as opposed to the theoretically favored 'bottom up' scenario). 

Because the power spectrum of density fluctuations in the Universe at the time of 
recombination is tilted to low masses in the ACDM scenario, hierarchical accretion is a 
necessary consequence of the hypothesis that galaxies are formed within cold dark matter 
haloes. The evidence presented here poses a severe challenge to the hierarchical formation 
scenario in that the observations show the opposite behavior to the theoretical predictions, 
with the more massive galaxies being already present at a lookback time of 65% of the 
Hubble time. 

Note added in proof: D. Stern (priv. comm.) has obtained a deep Chandra image 
of QSO1215-00 and finds no evidence of diffuse X-ray emission. 
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Fig. 2. — Rest frame from our data (filled points) in the 3.6//m and 4.5//m bands (arbi- 
trarily shifted by 0.05 in z for clarity) , together with previous -fi '-band studies (as identified 
in the figure legend) and models from lBruzual fc CharlotI (120031 ) . with Zj qjs, indicated in the 
legend and r = . 1 Gyr . The models are normalized to the value of K* in the Coma cluster 



De Propris et all (jlQQSh . 



Table 1. Observed clusters 
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60 


Romer et al. 2000 


RDCS1634.5+5724 


16 


34 


27 


6 


+57 


22 


51 


8 





61 


Rosati et al. 1998 


RDCS0046.3+8531 


00 


46 


19 


7 


+85 


31 


01 


3 





62 




RDCS0440.5-1630 


04 


40 


28 


4 


-16 


30 


08 








62 


?; 


RDCS0542.8-4100 


05 


42 


50 


2 


-41 


00 


07 








64 




MS1610+66 


16 


10 


47 


8 


+66 


08 


41 








65 


Gioia et al. 1990 


RDCS1936.0-4640 


19 


36 


06 


6 


-46 


40 


03 


6 





65 


Rosati et al. 1998 


RDCS0047.3+8506 


00 


47 


14 


8 


+85 


06 


02 








66 


" 


RDCS2313.6+1415 


23 


13 


34 


5 


+14 


15 


15 


5 





67 


" 


RDCS2038.5-0125 


20 


38 


29 


1 


-01 


25 


11 


7 





68 


" 


RDCS2236.7-2609 


22 


36 


42 


7 


-26 


09 


30 








70 


" 


RDCS2303.7+0846 


23 


02 


47 


5 


+08 


44 


07 


4 





73 


" 


GHO1322+30 


13 


24 


48 


2 


+30 


11 


14 








75 


Gunn, Hocssel & Oke 1986 


RDCS1517.9+3127 


15 


17 


56 


3 


+31 


27 


27 








75 


Rosati et al. 1998 


MS1137+66 


11 


40 


22 


3 


+66 


08 


15 








78 


Gioia et al. 1990 


RDCS1350.8+6007 


13 


50 


46 


1 


+60 


07 


09 


5 





80 


Rosati et al. 1998 


RDCS0035.9+8513 


00 


35 


55 


2 


+85 


13 


20 








81 


55 


RDCS1317.4+2911 


13 


17 


21 


4 


+29 


11 


25 








81 


5? 


RJ1716+67 


17 


16 


49 


6 


+67 


08 


30 








81 


Henry et al. 1997 


RDCS0152. 7-1357 


01 


52 


43 


7 


-13 


57 


21 








83 


Rosati et al. 1998 


RDCS0337.4-3457 


03 


37 


24 


7 


-34 


57 


29 








84 


55 


RJ1226+33 


12 


26 


54 





+33 


32 


00 








89 


Ebeling et al. 2001 


GHO1604+4304 


16 


04 


23 


2 


+43 


04 


44 








90 


Gunn, Hoessel & Oke 1986 


3C184 


07 


39 


24 


5 


+70 


23 


11 





1 


00 


Dcltorn et al. 1997 


MG 2019.3+1127 


20 


19 


18 





+11 


27 


10 





1 


00 


Hattori et al. 1997 


RDCS0910.7+5422 


09 


10 


45 





+54 


22 


02 





1 


11 


Rosati et al. 1998 


3C210 


08 


58 


09 


9 


+27 


50 


52 





1 


16 


J.-M. Deltorn, priv. comm. 


3C324 


15 


49 


48 


9 


+21 


25 


38 





1 


21 


Dickinson 1995 


RDCS1252.9-2927 


12 


52 


54 


2 


-29 


27 


07 





1 


24 


Rosati et al. 1998 


RDCS0848.9+4452 


08 


48 


56 


2 


+44 


52 


00 





1 


26 




RDCS0848.6+4453 


08 


48 


34 


2 


+44 


53 


35 





1 


27 




QSO1215-00 


12 


15 


49 


8 


-00:34:34.0 


1 


31 


Liu et al. 2000 
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Table 2. Number counts in clusters 



Cluster Mi^ui ^total ^background ^stars ^cluster ^lim ^total ^background ^stars ^cluster 





3.6/im 3.6/im 


3.6/im 


3.6/im 


3.6/im 


4.5/im 4.5/im 


4.5/um 


4.5/im 


4.5/im 


RJ1120+43 


17, 


.34 


120 


34.6 


6.4 


80.0 


17, 


.32 


168 


71, 


,4 


6.4 


90.2 


CL1634+57° 


17, 


,36 


75 


48.0 


14.4 


12.6 


17, 


,37 


112 


74, 


,8 


22.5 


14.7 


CL0046+85° 


17, 


,37 


131 


52.4 


35.5 


43.1 


17, 


,40 


177 


94, 


,6 


39.0 


43.4 


CL0440-16 


17, 


.23 


115 


47.8 


15.9 


51.3 


17. 


,40 


138 


72, 


,0 


16.5 


49.5 


CL0542-41 


17, 


.40 


161 


47.3 


26.5 


87.2 


17, 


,43 


189 


74. 


,7 


22.3 


92.0 


MS1610+66 


17, 


,45 


121 


49.1 


12.3 


59.6 


17, 


,44 


170 


73. 


,2 


14.2 


82.6 


CL1936-46'' 


17, 


,45 


134 


53.0 


89.3 


-8.3 


17. 


,44 


164 


75, 


,1 


89.2 


-0.3 


CL0047+85" 


17, 


,49 


133 


53.0 


38.2 


44.2 


17, 


,53 


147 


75, 


,1 


29.3 


42.6 


CL2313+14 


17, 


,53 


123 


58.5 


14.6 


49.9 


17, 


,54 


148 


84, 


,8 


16.5 


46.8 


CL2038-01 


17, 


,57 


183 


58.9 


14.7 


109.4 


17, 


,58 


191 


87, 


,4 


16.6 


87.0 


CL2236-26 


17, 


,64 


121 


58.7 


10.4 


51.9 


17, 


,62 


150 


88, 


,4 


16.4 


50.0 


CL2303+08" 


17, 


,27 


73 


24.0 


12.8 


36.2 


17, 


,69 


126 


88, 


,0 


20.7 


25.3 


CL1322+30 


17, 


,83 


134 


74.0 


5.2 


54.7 


17, 


,71 


157 


92, 


,2 


2.2 


62.6 


CL1517+31" 


17, 


,67 


105 


59.0 


9.7 


36.3 


17, 


,71 


117 


69, 


,4 


9.9 


37.7 


MS1137+66 


17, 


,74 


126 


63.7 


7.1 


55.2 


17, 


,77 


167 


94, 


,0 


7.9 


65.1 


CL1350+60 


17, 


,78 


139 


63.7 


6.3 


69.0 


17, 


,79 


181 


95, 


,3 


7.0 


78.7 


CL0035+85 


17, 


,81 


154 


64.4 


28.2 


61.4 


17, 


,82 


177 


96, 


,3 


30.9 


49.8 


CL1317+29 


17, 


,81 


114 


64.7 


5.3 


44.0 


17, 


,82 


139 


95, 


,9 


5.2 


37.9 


RJ1716+67 


17, 


.81 


167 


64.4 


19.8 


82.8 


17. 


,82 


184 


97. 


,3 


22.2 


64.5 


CL0152-13 


17, 


.85 


161 


67.1 


6.0 


87.9 


17. 


,82 


189 


97. 


,5 


6.8 


84.7 


CL0337-34" 


17, 


,88 


106 


67.7 


12.8 


25.5 


17. 


,84 


144 


97, 


,0 


14.6 


32.4 


RJ1226+33 


17, 


.97 


129 


70.5 


4.7 


53.8 


17, 


,85 


143 


93, 


,5 


5.2 


44.3 


GHO1604+4304'' 


17, 


.97 


112 


69.9 


11.6 


30.5 


17, 


,91 


132 


96, 


,2 


11.8 


24.0 


3C184° 


17, 


,16 


72 


57.7 


11.2 


3.1 


17, 


,02 


51 


34, 


,6 


11.3 


5.1 


AXJ2019.3+1127" 


17, 


,16 


146 


69.8 


75.6 


0.6 


17, 


,02 


134 


35, 


,0 


98.0 


1.0 


CL0910+54 


17, 


,37 


59 


33.7 


5.8 


19.5 


17, 


,12 


72 


38, 


,7 


5.1 


28.2 


3C210 


17, 


,32 


62 


29.1 


20.7 


12.2 


17. 


,18 


74 


40. 


,7 


6.1 


27.2 


3C324 


17, 


,43 


71 


35.1 


8.5 


27.4 


17, 


,24 


73 


41, 


,6 


8.7 


22.7 


CL1252-29 


17, 


,41 


102 


35.1 


15.8 


39.0 


17, 


,27 


93 


45, 


,9 


16.0 


37.1 


CL0848.9+4452'' 


17, 


.54 


39 


8.6 


1.7 


34.7 


17. 


,30 


21 


11. 


,1 


1.7 


9.3 


CL0848.6+4453 


17, 


.54 


69 


34.4 


6.5 


28.1 


17. 


,30 


79 


46, 


,3 


6.7 


26 


QSO1215-00'' 


17, 


.75 


55 


46.2 


4.8 


4.0 


17. 


,36 


48 


45. 


,4 


4.4 


-1.8 



Low number counts 
^ Only 0.5 Mpc field - cluster is off centre 
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Table 3. A;— corrections 



Redshift 


^3.6 


^4.5 


0.01 


-0.216 


-0.216 


0.025 


-0.249 


-0.268 


0.05 


-0.3 


-0.352 


0.075 


-0.347 


-0.433 


0.1 


-0.391 


-0.512 


0.125 


-0.434 


-0.588 


0.15 


-0.475 


-0.66 


0.175 


-0.516 


-0.724 


0.2 


-0.555 


-0.779 


0.225 


-0.594 


-0.833 
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Table 4. Luminosity Function Parameters 



Redshift 


M* (3.6;um) 


a 


M* (4.5//m) 


a 


0.75 


-24.63 ± 0.22 


-0.25 


-24.48 ±0.31 


0.21 


1.15 


-26.18 ±1.20 


-0.82 


-24.83 ± 0.83 


0.81 



